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We shall present here the causal interpretation of canonical quan- 
tum gravity in terms of new variables. Then we shall apply it to the 
minisuperspace of cosmology. A vacuum solution of quantum cosmol- 
ogy is obtained and the Bohmian trajectory is investigated. At the 
end a coherent state with matter is considered in the cosmological 
model. 
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1 Introduction 

The causal interpretation, introduced by de-Broglie|lj and Bohm 2J, presents 
a definite trajectory for any quantum system. It is proved tliat Copenliagen 
and causal interpretations answers to any common question are the same. 
There are questions specific to causal interpretation, directly related to the 
trajectories which can not be asked in Copenhagen interpretation. Thus it 
is not simply another way of interpreting the same theory. But the problem 
that whether such questions are physical or not is still unclear jS]. 

The causal interpretation is motivated by the fact that the phase of 
the wave function for a non relativistic particle obeys the Hamilton-Jacobi 
equation modified by a potential called quantum potential. Interpreting this 
equation as a quantum Hamilton-Jacobi equation one can derive the quan- 
tum trajectory. The power of this interpretation comes to the scene when 
one considers problems like measurement. It gives a deterministic version 
of measurement theory with the same results as the Copenhagen interpre- 
tation. Extension of the theory to the case of relativistic particles jl] and 
fieldsjS] is straightforward. It is also possible to make a causal interpretation 
of quantum gravity in terms of old variables jS]. 

Here we shall investigate the possibility of making a causal interpreta- 
tion of quantum gravity in terms of new variables [7J (loop quantum gravity). 
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Causal loop quantum gravity may be useful at least in three points. The 
first point is the issue of time. In loop quantum cosmology the notion of 
time can be introduced either internally, or externally. In the first view the 
scalar constraint is looked at as an evolution equation with respect to the 
eigenvalues of momentumjHj. In the second view one introduces an external 
time which parameterizes quantum theory[9 . It has no physical meaning 
but provides time dependent states, called state-time^. The time issue in 
Bohm's theory is different from these two definitions of time. In Bohmian 
mechanics the guidance relation connects canonical momentum to the gradi- 
ent of the Hamilton-Jacobi function. Therefore one has a definite trajectory 
parameterized by this time, which is exactly the physical time for classical 
trajectory. In this case there is no approximation and if matter exists, one 
can eliminate the coordinate time, using the evolution of the matter fields, 
and introduce a clock field. 

Another point is the evolution of the universe through the big bang singu- 
larity. It is shown that in loop quantum cosmology the quantum evolution 
is well defined through the big bang singularity. Investigation of existence or 
non-existence of singularity in causal loop quantum cosmology, means look- 
ing for any possible singularity in Bohmian trajectory. It is a general feature 
of de-Broglie-Bohm theory that it is quiet possible to have singularity in the 
wave function without any singularity in the trajectory. 
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The last point is about the role of the wave function. Causal quantum 
theory can talk about single systems, like universe, without any reference to 
the probabilistic character of the wave function. Also in this theory, there 
is no need to have a classical system outside the observed system (the ob- 
server) to interpret the measurement problem. This shows the usability of 
the application of de-Broglie-Bohm theory to quantum cosmology. 

In the next section we drive the guidance relation for loop quantum grav- 
ity and in section 3 we shall evaluate it for the cosmology minisuperspace. 
In section 4 we shall find exact solutions and their Bohmian trajectories. 

2 Causal loop quantum gravity 

Loop quantum gravity is based on the formulation of general relativity in 
terms of su{2) connections (^4^) and triads {Ef). Einstein's equations are 
expressed in terms of three constraints. Two of them (gauge and 3-diff 
constraints) are kinematical and the last one (the scalar or Hamiltonian con- 
straint) is dynamical. 

In loop quantum gravity the elementary variables are constructed from 
the above fields. They are holonomies along edges of a graph and smeared 
fluxes of triads through surfaces. These have quantum representations in 
the Hilbert space of the functions of holonomies. An orthonormal basis for 
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this kinematical Hilbert space is spin networks which are eigenvectors of 
geometric operators and depends on the connection via holonomies. 

In order to have a causal interpretation of any quantum theory, one needs 
a guidance relation as well as the wave function. The guidance relation en- 
ables one to derive Bohmian trajectories from the phase of the wave func- 
tion. For the case of loop quantum gravity the wave functional is obtained 
by solving the above mentioned three constraints. The guidance relation can 
be obtained by using the Poisson bracket between holonomies and fluxes. 
Therefore if we denote the holonomy along an edge e as A{e) and the elec- 
tric flux through a surface A as P( A, /) = /jS' (S^^ = rjabcP^^ is the 
2-form dual of the electric field = E^/STcG'-f, 7 is the Barbero-Immirizi 
parameter and / is a test function.), we have the commutation relation |7j: 



where k, equals to +1(-1) if e lies above (below) A and zero otherwise and 
p is the point that the surface A intersects with the edge e. Let us now 
derive the guidance equation. First consider a simple case when we have 
a Poisson bracket like {q,p} = X{q). The quantum system would have 
the commutation relation [q,p\ = ihX{q), leading to the q-representation 
p = -ihXd/dq and thus the Bohmian relation 

PBohm 

XdS/dq. The 
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extension to our case would be: 

trace (^^j^y4(e)r*/i(p)) if P is the source of e 
trace (-^|^/i(p)rM(e)) if p is the target of e 

(2) 

Thus in order to derive the Bohmian trajectories one should first solve the 
three quantum constraints and obtain the wave functional. Then denoting h 
times the phase of the wave functional as S, we can use the above guidance 
relation to obtain the quantum dynamics of the geometry. We shall do this 
for cosmological minisuperspace in the following section. 



3 Causal loop quantum cosmology 

Using the spatial homogeneity and isotropy of cosmological models one can 
find solutions of the 3-diff and Gauss constraints [S]. The elementary variables 
are holonomies along straight edges and fluxes across squares. The kinemati- 
cal Hilbert space is made from gauge invariant isotropic states which depends 
on the connection via the holonomies and form an orthonormal basis. 

In such a minisuperspace the connection has only one degree of freedom 
c = ^{k — 7d) where a is the scale factor and k = 0,+l is the curvature 
parameter for a flat and closed model respectively, and the conjugate mo- 
mentum is p with \p\ = a?. 
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The action of scalar constraint after regularization and quantization on 
physical states leads to the following difference equation Ej: 

where VL = 2 — 4£Q7^r(r — 1), £ is the eigenvalue of p, 4 is a quantum 
ambiguity coming from regularization |7^, Cmatter is the matter Hamiltonian, 
and ip is the Planck length. T is the spin connection related to the curvature 
parameter such that F = 0, 1/2 for the flat and close universe respectively. 
ip{4>, i) are coefficients of the expansion of the state in terms of spin network 
states: 

\iJ) = Y^^P{<P,i)\i) (4) 

e 

in which represents matter fields and = ^ £3 jg ^.j^g eigenvalue 

of the volume operator. In order to the sum make sense, one should add the 
condition that i) has support only on a countable subset of real numbers. 
The gauge invariant isotropic state can be represented by {c\i) = e^^^ which 
form an orthonormal basis. The presence of Iq is related to the quantum 
ambiguity and it can be fixed by minimal area of the full theory. 

As one may expect, the discrete evolution in terms of internal time ^ 
does not break at big bang, this leads to resolution of singularity problem 
in quantum cosmology[8j. At scales enough larger than 4? the differential 
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WDW equation emerges as an approximation 3,- Also there is a lot of works 
done on effective field equation which incorporate quantum corrections to 
FRW equations [Ej. 

In order to have a causal interpretation, one needs the guidance relation. 
It can be derived from the guidance relation of the full theory given by the 
relation (0) by going to the minisuperspace of cosmology. Taking /j = r^, 
since the holonomy for an edge of length I is A{e) = cos ic/ 2 + 2r* sin^c/2, 
using the relation (j21), we have: 

Tc=P 

In the next section we shall solve for the wave function in the absence of 
matter and derive the Bohmian trajectory. 



4 Cosmological solutions 
4.1 Vacuum solution 

For simplicity here we shall restrict ourselves to vacuum solutions of equation 
(jHI). In the vacuum case the difference equation (jH} can be written as: 

F{i + 4io) - nF{i) + F{e - Mo) = (6) 

where F{i) = {Ve+e,, - Vg^e,,)^{i)e-'^^/^ . Note that for £ = this leads to 
F = and thus does not determine iplO). Using the Z-transform techniques, 
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considering a solution of type F{i) ~ e^^, the above equation leads to: 



P={ " • (7) 

Since the Z-transform solution in the case /c = is degenerate, one has a 
second solution of the form ie^^. Therefore the general solution of equation 
© is: 

3/2^ I e^^''^^/2^o(A + 5£); k = 



m 
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Ae*^ + 5e-^o^; A; = 1 
where A and -B are two constants. Note that this is only for i ^ 0, and that 
?/'(0) is just another free parameter. 

For large i, we have the following asymptotic behavior: 



~ J\e\; fork = 0,A = (},B ^0 (9) 



^(i) ^e^o\e\. fork = l 



The behavior of the wave function is plotted in figures (0)- Q. It is clear 
from these figures that only for the case k = and B = 0, the wave function 
converges for large i. The divergence in the other cases, is quiet similar to the 
results of jTTj. Here, thus we choose the convergent case. As £ is a continuous 
variable, one can construct the general solution as the following form: 

^ijTr(n+no)/2£o 

noG[0,l) neZ 
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where /(no) is some coefficient, and again has support only on a countable 
subset. In the above relation we have broken the sum over the real number 
£ into sum over a unit interval and an integer n. The term no = in the 
above summation has the problem that setting n — leads to an infinite 
contribution. This suggests that one needs to introduce a condition on /(no) 
as /(O) = and this makes a finite state. 




Figure 1: The real (solid line) and imaginary (dashed line) parts of the wave 
function as a function of i. This is plotted for the case k = 0, A = 1 and B = 0. 
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Figure 2: The real (solid line) and imaginary (dashed line) parts of the wave 
function as a function of i. This is plotted for the case k = 0, ^ = and B = 1. 

4.2 Bohmian trajectories 

In order to get some physical details from the wave function ()10|). we restrict 
ourselves to the case that f{no) is equal to zero for all no except a specific 
one, no- To derive the Bohmian trajectories, we use the wave function 
in the configuration space. One can get an approximate closed form for the 
physical wave function noting that for large \n\ (i.e. grater than |no + ^o|): 

\n + no + iof^^ - |n + no - ^o|^^^ = 3£osgn{n)^J\^\ + (11) 
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Figure 3: The real (solid line) and imaginary (dashed line) parts of the wave 
function as a function of i. This is plotted for the case k = 1, A = 1 and B = 0. 

Introducing <^ = —jj^^ — c, the wave function is approximately: 

{-| -| oo —ina i oo ina 
- + Lyi I 

(12) 

The last two terms can be summed to an integral form|12j: 

i sin a r°° dx i 



/■oo {i-r 9 r ■ ■ ■ ■ 

/ = -- <|.(e-", 1/2, 1) - e^-$(e*", 1/2, 1) 

^ Jo ya;(coshx — cosaj 2 l 

(13) 

where $ is the Lerch (-function ^{z,s,v) = (T{s))-^ dt{f-^e-''^) / {I - 
ze~^). 
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Figure 4: The real (solid line) and imaginary (dashed line) parts of the wave 
function as a function of i. This is plotted for the case k = 1, ^ = and B = 1. 

To drive the Bohmian trajectory one needs to extract the phase of the 

wave function in the configuration space. This can be done numerically. 

The result is plotted in figure (0). In figure © the Bohmian trajectory 

for the scale factor is plotted. In these figures we have used io = VSn, 

7 = 0.238 and j = 0. As it is seen in figure (0), the quantum Hamilton- 

Jacobi function has some oscillations around the classical Hamilton-Jacobi 

function {Sc = —hoa)} As a result, Bohmian trajectory in figure ^ is also 
^This can be simply understood. Since the classical solution of such a model (FRW 

without matter) is a static universe (a ^constant), we have dSc/Sc ^constant and thus Sc 
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oscillatory around the classical path, which is a constant scale factor in this 
case. 




Figure 5: Hamilton- J acohi function (solid line is the quantum case and dashed 
line is the classical case) plotted versus a. 

4.3 Inclusion of matter 

For simplicity we assume that the matter Hamiltonian is such that, an initial 

state which is highly peaked at the classical solution, would be remained so 

for a long time. This condition is satisfied by choosing a Gaussian distribution 
is a linear function in c. 
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Figure 6: Classical and quantum trajectories of the scale factor. 

like: 

m = Y.HnA)e- = exp (-^^) (14) 

Here c is the classical connection and a is a complex function determining 
the width of the Gaussian, and we denote its real and imaginary parts as: 

(15) 



a^{ct>) al{ct>) al{<p) 

The scalar constraint (jHl), determines the dependence of width a of the Gaus- 
sian distribution on matter. Here we shall assume that the state is so peaked 
at the classical path that one can make a classical approximation for matter 
field (setting (p = (p where (p is the classical matter field.). That is to say one 
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can substitute matter fields in the guidance relation of the scale factor with 
its classical value. The Hamilton- Jacobi function is then given by:^ 

^ = -^ (16) 

Using the guidance relation and the fact that \p\ = a?, the quantum path of 
the scale factor is given by the equation: 

a = S±^a^ (17) 

7 

Therefore for such a semi-classical state the Bohmian trajectories are fluc- 
tuations around the classical path. 
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